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Abstract 



We study the effect of Aharonov-Bohm flux on tlie superconducting state 
in metallic cylinders. Although Byers and Yang attributed flux quantization 
to the flux-dependent minimum of kinetic energies of the Cooper pairs, it is 
shown that kinetic energies do not produce any discernible oscillations in the 
free energy of the superconducting state (relative to that of normal state) as 
a function of the flux. This result is indeed anticipated by the observation 
of persistent current in normal metal rings at low temperature. Instead, 
we have found that pairing interaction depends on the flux, leading to flux 
quantization. When the flux ($) is given by $ = n x hc/2e (with integer 
n), the pairing interaction and the free energy become unchanged (even n) 
or almost unchanged (odd n), due to degenerate-state pairing resulting from 
the energy level crossing. As a result, flux quantization and Little-Parks 
oscillations follow. 



PACS numbers: 74.20.Fg, 74.25. Jb, 74.62.-c 



1 



I. INTRODUCTION 



Flux quantization is one of the most fascinating properties of superconducting states. 
This phenomenon was found by Deaver and Fairbank^ and Doll and Nabauer^ in 1961. They 
observed that the magnetic flux trapped in a superconducting cylinder is quantized in units 
of hc/2e. Recently, half-integer flux quantization has attracted much attention in connection 
with the pairing symmetry in high Tc cuprates.'^ 

Byers and Yang^ proposed that flux quantization follows because the free energy of the 
superconducting state is periodic in the unit of the flux, hc/2e, if electrons are paired. In 
fact, they assumed that i) the free energy, F^, of the normal state is independent of the 
flux, i.e., Fn{^) = constant, and ii) the pairing energy is constant as a function of the flux.^ 
Consequently, they ascribed flux quantization to the minimum of the kinetic energy of the 
superconducting state at $/(/ic/2e) = integer. Based on the theory of Byers and Yang, 
Little and Parks^ pointed out that the transition temperature Tc is also a periodic function 
of the enclosed flux which is called Little-Parks oscillations. 

However, recent discovery of persistent current in normal metal rings at low temperature 
casts a doubt on the first assumption, that the free energy of the normal state is essentially 
independent of the flux.^~^ Since the electronic states of normal metal rings are periodic 
in flux with period hc/e, the persistent current is given by /($) = — c x 9F7v($)/9$. 
Therefore, the first assumption is not compatible with the experimental observation of the 
persistent current in normal metal rings. 

In this paper we report that fiux quantization and the Little- Parks oscillations are caused 
by the fiux dependence of the pairing energy. In addition, we show that the fiux dependence 
of the kinetic energy does not lead to any significant oscillations in the transition temperature 
Tc and the energy gap parameter Aq for superconducting cylinders. As the fiux increases, 
the pairing matrix element is reduced since we must pair single-particle states which are 
not degenerate and have different density distributions. ^"'^^ When the fiux ($) is equal 
to $ = n X hc/2e (with integer n), the pairing interaction and the free energy of the 
superconductor remain unchanged (even n) or almost unchanged (odd n) from their values 
with zero fiux, since we can pair degenerate partners, with the same density distribution, 
due to the energy level crossing. Accordingly, we find fiux quantization and Little-Parks 
oscillations. It is interesting to note that Schwartz and Cooper^^ obtained a similar result 
for generalized pairing with different angular momentum in superconducting cylinders. 



II. THEORY OF BYERS AND YANG 

To study the effect of the Aharonov-Bohm fiux $ on the superconducting states in 
cyhnders, Byers and Yang^ considered the variation of the energy levels E of a single particle 
in the presence of the fiux: 

Here Pr and p^ denote momenta in the radial and z directions. Employing the BCS pairing 
between states m and -m, they noted that the average (kinetic) energy for the states m and 
-m, near $ = -|-0, increases with $ like 
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constant + - — ^ t- • 2) 
2mr^ he 

Furthermore, at 2e^/hc — 1, sinee pairing between 

e ^ 11,13,3 
m+— $ = m + - = - and--, - and--, etc., (3) 

is preferred, the kinetic energy per particle remains the same as for the case $ = 0. Then 
the additional energy for each of these pairs, in the neighborhood of 2e$//ic, increases with 
$ in a parabolic manner: 

, xr/C^x 1,9 

Consequently, the kinetic energy of the Cooper pairs has the minimum at $ = n x 
hc/2e (with integer n). Based upon this observation, Byers and Yang^ asserted that flux 
quantization is due to the flux-dcpcndcnt minimum of the kinetic energy of the Cooper pairs. 
On the other hand, the corresponding BCS wavefunctions at $ = n x hc/2e are given by^^ 

vOi^-^-mi^^v+mi^l^ ^1 even) (5) 

n(iii.+i/2 + v^+ii2a^^^^^^at^^^^)\Q >■ {n odd) (6) 

with n — 2m {n — 2m + 1) for even (odd) values of n. The operator a^^_^| creates an 
electron in the state {v -\-m\). 

However, the kinetic energy of the normal electrons will show the same flux dependence, 
as confirmed by the observation of the persistent current in normal metal rings. Therefore, 
we expect that this fiux dependence will be canceled in the superconducting state. More 
precisely, using the BCS theory, we show that the variation of the kinetic energy with the 
flux does not lead to any oscillations of the condensation energy (i.e., Aq) and the free energy 
(i.e., Tc) of superconducting cylinders. 

Assuming a constant pairing matrix element V,^'^ one finds the BCS gap equation at 
T = Ois: and the equation:^^-^^ 

2^( '-^ )2 + A^ 



m' 



tanh^^^^^ + tank- 



where e„i-^ is the one-particle energy (measured from the Fermi energy Ep) with = e^/hc. 
Contrary to Byers and Yang, it is clear that the kinetic energy does not give rise to any 
significant oscillations in the gap parameter, A^ = Aq, since its effect is only to shift the 

t; 2 2 

BCS cutoff region from —ujd < e < ud to —ujd + < e < (^d + (where = 2~^)- 



Ap^-yy , ^° (9) 

m 2^ {em + e^Y + Ag 
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Use of em-n + ^-m-fj, = 2(e,n, + e^) has been made. For the Tc, Eq. (8) shows indeed a 
minor decrease in Tc (for nanoscale systems) due to the flux, which is inconsistent with 
the (T = OK) gap equation, Eq. (7), and has been elusive.^'® However, this decrease is 
completely neghgible for bulk systems and it originates from the exclusion of the Cooper 
pairs with negative excitation energy. ^^'^^ Since analytic solution is not available, we rely on 
the numerical calculations and merely emphasize that the correction is in the denominator: 

1 = -F ^ ^ 27^ 

^ 1 sinh'-^ ^ ^ 

^-VY^ ^ e_^^ (10) 

^(^772 ~l~ ^/j-} cosh 2^ cosh 2^ 

Note that for bulk systems m is of the order of 1,000 near Ep, while /x is of the order of 1. 

For example, we solve Eqs. (7) and (8) numerically for a superconducting cylinder with 
the inner radius Va — 7, 200A, the outer radius = 8, OOOA, and the height L — 500A. 
For a thin cylinder, i.e. ri, — Va << H, the electron energy is approximately written as 
em-M = £[^^7177^ + iv^y + (^)2].^9'20 (More calculational details will be explained 
in the next section.) We assumed Ep = l.OeV and up, = 0.026^. There are 1,286,200 
states in the BCS cutoff range. Figure 1 shows Aq and versus n = e^/hc. As can be 
seen, there is no change in Aq (i.e., the condensation energy) and Tc (i.e., the free energy 
of superconducting state) due to the variation of the kinetic energy with the flux For a 
certain V (= 17.0/volume [eV]), the gap parameters Aq are the same 2.144 886 meV for = 
and /i = 0.25, whereas the T^s are 1.212 888 meV {fi = 0) and 1.211 921 meV (/i = 0.25), 
respectively. It is noteworthy that the decrease in Tc scales approximately ~ 1/ ^/N, with 
N= number of states in the BCS cutoff range. Consequently, the Tc decrease due to the 
exclusion of the Cooper pairs with negative excitation energy^^'^^ can be completely ignored 
in the bulk limit. Note also that 2Aq/Tc = 3.537, in agreement with the bulk result. 



III. FLUX-DEPENDENT PAIRING INTERACTION 

Now we consider flux-dependence of the pairing matrix element and the pairing interac- 
tion in a superconducting cylinder whose axis coincides with the z-axis. The inner and outer 
radii are and r^, and the height is L. First, we determine one-particle eigenfunctions and 
eigenvalues in the presence of the Aharonov-Bohm flux The Hamiltonian is given by 

where /j, = e^/hc. Assuming a solution of the form 

0(r, 0, z) = R{r)sinC^)e''^'f' . (12) 
we obtain the Bessel equation for R{r): 

R" + ^ + [k' - ^"""/^V = 0. (13) 
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Prom the boundary condition R{r — Tq) = R{r — r^) = 0, we find the radial solutions and 
eigenvalues:^^'^° 

Rm-A^) = ^^[Nrn-^,{kra)Jm-^^{kr) - J,n-,.{kra)Nm-^{kr)] (14) 

Eimn. = + (^)'] (n, = 1, 2, 3, • • • , m = 0, ±1, ±2, • • •) (15) 

where H is the normalization constant and kmi is the Zth root of the transcendental equation 
N^_^{kra)Jm-n{kn) - Jm-n{kra)N^_^,{kn) = 0. (16) 

Prom Eq. (16) wc calculate the single-particle energy as a function of the fiux shown 
in Figure 2, near Ep = 1.0 eV for = 75 A and Vb = 150 A. Ignoring the z coordinate, 
we find 15 basis pair states within the BCS cutoff range for ujd = 0.02eV. We consider 
a small size ring for simplicity and because of the difficulty in calculating the high order 
Bessel functions. As /j, increases, some states near 1.02 eV get out of the BCS cutoff range, 
while some states below 0.98 eV are coming in the range. We disregard this complication, 
since this effect will be negligible in the bulk limit. Moreover, it is caused by the simple 
approximation, which can be avoided in the strong-coupling theory.^^ Near n = 0.5 we chose 
Ep = 0.99 eV, which is obtained from the Ep for /i, = by applying the flux. Although the 
spectra near /i — and /i — 0.5 are similar, they are not exactly the same due to the finite 
size effect. Notice the crossings of energy levels at ^ — hc/2e x n (with integer n). 

Second, we derive the pairing matrix element for an Einstein phonon model with the 
phonon Green's function D[x — x'). Using the equivalent electron-electron potential in the 
electron-phonon problem,^^'^^'^^ 

V{x-x')^VqD{x-x'), (17) 

with a; = (r, t), we get the pairing matrix element: 

Vm,m' = Vo 12 Rl,^^{r)R*_^,_^{r)R_^_,{r)Rrn^^{r)rdr; -1 + 5 < < ^ - 5 (18) 

Vm,m' = Vo r Rl,,_^{T)R*_^,^^_^{r)R_^^,_,{T)Rrn-,{T)TdT- ]- - 5 < I, < I - 5 (19) 
Jra ^ 

rn 3 
Vm,m' = Vo / R*^,^^ {r)R*_^,^^ {r)R_m+i-ix{r)Rm+i-ix{r)rdr] 1-5 <n<--5 (20) 

where 5 ~ 1/4. We have made use of the BCS wavefunction with the flux: 

(t>BCS = HiUm + ^^ma;^-^T«-m-MjO ^ + ^ < < \ - ^ (21) 



(pBcs = Y[{um + Vma^_|^^atrn+l-^,i\'^ - -S<fi<l-5 (22) 

m ^ 
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^Bcs = '[[{um + 'yma;^-^tal,rz+2-/.i |0 >; l-S<n<--5. (23) 



Finally, we obtain the BCS gap equation 
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where the excitation energy Em-ij. is 



Em-, = lie^-, - e^m-,) + + Al. (25) 

Since the analytic expression of the matrix element is not available, we rely on the numerical 
method to calculate from Eq. (24). Note that Eq. (24) leads to the negative excitation 
energy, when e^-ju and e-m-n have different signs near Ep, as in the cases of a supercon- 
ductor with a uniform exchange field^^ and the current-carrying superconducting state. 
Nevertheless, in this case the possibility of the negative excitation energy is not significant 
as shown in Fig. 1, because m ~ 1,000 near Ep while ~ 1 in the bulk limit. In other 
words, the exclusion of the contribution of the negative excitation does not lead to any 
noticeable effect in T^, as shown in Fig. 1. 



IV. EXPLANATION OF FLUX QUANTIZATION 



Flux quantization may be explained from the flux-dependence of the pairing matrix 
element Vmm' ■ We stress the periodicity of the matrix element with the flux: 

he 

Vmm'i^ = ^ X even n) = V^^, ($ = 0) . (26) 



V^mi^ = ^xoddn)=V^^,{^ = Q). (27) 

This behavior is caused by the energy level crossings at $ = hc/2e x n {or fj, — 
integer or half — integer). It was shown that the level crossings (at integer /i) are due to the 
cylindrical symmetry of the system.^'^'' Since for odd n the energy level crossing occurs at 
lower energy than for even n, the matrix elements are slightly different for even and odd n's. 

In the neighborhood of $ = 0-|- we pair states (m — /x |) and (— m — i) which are 
not degenerate and have a different density distribution. Accordingly, the pairing matrix 
element Vmm' is decreasing with the flux, compared to Vmm' at zero flux. For a large cylinder, 
we may employ the asymptotic radial wavefunction:^'^ 

Rm-,{r) = + am-,— - bm-A^, + ^)] ^^^} (28) 

where 
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_ 4(m-/i)^-l , _ [4(,»_^)2_i][4(^_^)2_9] 

and 

= (30) 
r-ft - r„ 

It is obvious that the algebraic terms lead to the different density distributions of states 
(m — /i t) and (— m — i) and the decrease of the matrix element: 

K.^,($) [1 - C"(^)V][l - C"(^)V] X yw($ = 0), (31) 

where C ^ (l/r^r;, - l/2rf - l/2rl)\ 

Figure 3 shows the radial density distributions |i?(r)pr of the m = ±54 states when 
/X = 0.25, (i.e., m — fi = 53.75, —m — /i = —54.25) for the same ring as Fig. 2. The solid line 
corresponds to m — /x = 54.75, whereas the dotted line corresponds to —m — = —54.25, in 
which the electron is rotating faster, due to bigger angular momentum projection along the 
z-axis, leading to the shift of the radial distribution towards the outer radius. The difference 
in the radial density distributions then induces the decrease of the pairing matrix element 
with increasing ^. Figure 4 shows the matrix element Vmm' = Vij for /i = 0. Notice that it 
is 15 X 15 symmetric matrix, since there are 15 pair states in the BCS cutoff range. The 
decrease of the matrix element with /i is not big enough to be noticeable in Fig. 4, though. 
We just mention the decrease of biggest matrix element Vm=53,m'=5i = l^=2j=3, between pair 
states with the same / value, from 0.0002450 to 0.0002448. 

This observation and the periodicity of Vmm' = Vij prove that flux quantization is due to 
the minimum of the pairing energy of the superconducting cylinders as a function of the flux. 
Accordingly, the free energy of the superconducting state also has minimum at $ = hc/2exn 
since the contribution of the kinetic energy is constant. In other words, pairing between the 
degenerate states caused by the energy level crossing at $ = hc/2e x n is the physical origin 
of flux quantization. Thus the Little-Parks oscillations follow accordingly. From Eqs. (24) 
and (31) we may estimate roughly the decreases of the energy gap at T = OK and T^. 

A^($) ~ A^($ = 0)(1 - 8C" < > ii'/NoVo), (32) 

T, ~ reo(l - 8C" < > //ViVoVb) , (33) 

where Nq is the density of states at Ep and <> denotes the average over the states in the 
BCS cutoff range. T^q is the transition temperature in the absence of the magnetic flux. 

Figure 5 shows the T = OK anisotropic gap parameter Aj for fJ. = 0. calculated from Eq. 
(24) by employing the matrix clement Vij = Vmm' as in Fig. 4. Even though the coupling 
between the second and third states is strongest, we find the biggest gap at the tenth state, 
i.e., Ai=io, due to its strong couplings to several neighbor states. 

Figure 6 displays the decrease of the (average) energy gap with yU. The solid line denotes 
the (T = OK) average gap parameter, A = J2i Ai/15. It should be noticed that this decrease. 
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so much enhanced as a result of the finite size effect, persists even in the bulk limit, as Eq. 
(32) shows. The different variations near ^ = and = 0.5 are also due to the finite size 
effect. The dashed line shows the isotropic gap parameter, according to the assumptions of 
Byers and Yang."^ It exhibits a small decrease with the fiux too, which quickly disappears 
with the increase of the system size. For the isotropic solution, we included Va — V besides 
Vij = V = constant. Therefore, we get the bigger gap. The semicircle is the Tc when fi = 0. 
Since, as the fiux is increased, the finite temperature gap equation, Eq. (24) eliminates the 
contributions from the Cooper pairs with negative excitation energy,^'''^^ drops to zero 
immediately. If we consider the Cooper problem, it is easy to understand why we need 
to eliminate the Cooper pair basis states with different signs in energy near Ep. However, 
this behavior is greatly amplified because of the finite size effect and it will not cause any 
problem in realistic systems, as shown in Fig. 1. Furthermore, the strong coupling theory^^ 
will reduce the effect even more. For nanoscale systems it may be desirable to use the 
canonical ensemble approach. We expect that for a reasonably large system, we find a 
continuous (small) variation of the Tc, consistent with the variation of the energy gap. 

V. DISCUSSION 

It is clear that more study is needed to understand the details of fiux quantization. In 
particular, it is essential to consider a big hollow cylinder for comparison with the experi- 
ment, which requires high order Bessel functions. This study indicates that the parabolic 
background^'^ in the Tc oscillations may be due to the magnetic field penetration into the 
superconductor, which is not addressed here. It is also interesting that disordered Auo ^Ino ^ 
cylinders^'' show half-integer fiux quantization as in the case of high cuprates.^ 

This study also implies that the exact calculation of the eigenenergies using high or- 
der Bessel functions and the magnetic field penetration into the metals are important in 
calculating the persistent current in normal metals. 

Finally, fiux quantization plays an important role in the vortex state^^ which may be 
understood by generalizing this approach. 

VI. CONCLUSION 

We have shown that fiux quantization and the Little-Parks oscillations in supercon- 
ducting cylinders are due to the fiux dependence of the pairing interaction, while the fiux 
dependence of the kinetic energy is canceled as evidenced by the observation of the persis- 
tent current in normal metals. When the flux $ is given by $ = hc/2e x n (with integer 
n), the pairing interaction and the free energy are unchanged (even n) or almost unchanged 
(odd n), due to the degenerate-state pairing, compared to those in the absence of the fiux. 
Accordingly, fiux quantization and the Little-Parks oscillations are obtained. 
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FIGURES 



FIG. 1. Energy gap and Tc as a function of the flux $ for a cylinder with Va = 7, 200A, 
Tf, = 8, OOOA, and L = 500A, following Byers and Yang. Energy gap does not change at all, while 
Tc shows negligible oscillations. 

FIG. 2. Onc-particlc energy levels vs ^ = e^/hc near the Fermi energy (Ep = l.OeV for /i = 
and Ep — 0.99eV for fi = 0.5) for a ring with ra = 75A and ri, = IbOA. Note the level crossings at 
fi = and 0.5. 

FIG. 3. Radial probability density \R{r)\'^r of states m — jjL = 53.75 (solid line) and 
— m — /X = —54.25 (dotted line) for the same ring as Fig. 2. Notice the shift of the radial 
distribution of —m — /J. = —54.25 state to the outer radius. 

FIG. 4. Pairing matrix element Vm^i = Vij for the same ring as Figs. 2 and 3. The biggest 
element Vi=2j=3 = Kn=53,m'=5i corresponds to the pair states with the same I value. The diagonal 
term Vmm = Vii does not contribute to superconductivity. 

FIG. 5. Anisotropic gap parameter Aj for = at T = OK, corresponding to the matrix 
element, Vij, as in Fig. 4. 

FIG. 6. Average energy gap parameter (solid line), A = X]i^i/15, (at T = OK) vs /x for the 
same ring as Figs. 2-5. The dashed line shows the isotropic gap parameter vs /x following Byers and 
Yang. The semicircle denotes the Tc for fi = 0, corresponding to the anisotropic gap parameter, 
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